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Abstract. We investigate the effects of a charge asymmetry on the spectrum of dileptons radiating 
from a quark gluon plasma. We demonstrate the existence of a new set of processes in this regime. 
The dilepton production rate from the corresponding diagrams is shown to be as important as that 
obtained from the Born-term quark-antiquark annihilation. 



INTRODUCTION 

The aim of this talk is to show that, when in a medium there is a finite charge density 
(i.e., a finite chemical potential), a new set of lepton pair-producing processes actually 
arises []. We then calculate a new contribution to the 3 -loop photon self-energy. The 
various cuts of this self-energy contain higher loop contributions to the usual processes 
of qq e+e^, qg qe^e^ , qq qqe^e^ , and an entirely new channel: gg e+e^. We 
calculate the contribution of this new reaction to the differential production rate of back- 
to-back dileptons. It is finally shown that, within reasonable values of parameters, this 
process may become larger than the differential rate from the bare tree level qq — > 

At zero temperature, and at finite temperature and zero charge density (note: hence- 
forth, a finite density will imply a finite charge density), diagrams in QED that contain 
a fermion loop with an odd number of photon vertices (e.g. Fig. P are cancelled by 
an equal and opposite contribution coming from the same diagram with fermion lines 
running in the opposite direction (Furry's theorem ^, ^]). This statement can also 
be generalized to QCD for processes with two gluons and an odd number of photon 
vertices. 

A physical perspective is obtained by noting that all these diagrams are are encoun- 
tered in the perturbative evaluation of Green's functions with an odd number of gauge 
field operators. At zero (finite) temperature, in the well defined case of QED we observe 
quantities like (0|A^iA^,2...A^,2n+i|0) ( Tr[p{fj,^)A^iA^2---A^2n+i] ) under the action of 
the charge conjugation operator C. In QED we know that CA^C^^ = — A^. In the case of 
the vacuum |0), we note that C|0) = |0), as the vacuum is uncharged. As a result 

(0|A,,A^,...A^,„^J0) = {0\C-'CA,^C-'CA,,...A,,„^^C-'C\0) 
= {0\A,,A,,...A,,„J0){-If"+' = -{0\A,,A,,...A,,„J0)=0. (1) 



Note that a net charge density in a quark gluon plasma does not necessarily imply a net baryon density 
and vice-versa. 



FIGURE 1. Diagrams rendered zero by Furry's theorem. See text for details. 



At a temperature T , the corresponding quantity to consider is 

n 

where P = l/T and is a chemical potential. Here, however, C\n) = e^^ \ —n), where 
I — n) is a state in the ensemble with the same number of antiparticles as there are 
particles in \n) and vice-versa. If ^ = i.e., the ensemble average displays zero density 
then inserting the operator C^C as before, we get 

{n\A^^A^,...A^,„Jn)e-^^" = -{-n\A^^A^,...A^,„J-n)e-^''\ (2) 

The sum over all states will contain the mirror term (— nlA^jA^j-'-^/iin+i I —n)e^^^", with 
the same thermal weight 

^ l^{n\A,^A,,...A,,„^^ = 0, (3) 

n 

and Furry's theorem still holds. However, if 7^ unequal number of particles and 
antiparticles ) then 

{n\A,,A,,...A,,„Jn)e-^(^"-^Q"^ = -{-n\A,,A,,...A,,„J -n)e-^^^"-^Q''\ (4) 

the mirror term this time is (— n|A^jA^2---^/i2n+i I — , with a different 

thermal weight, thus 

Y,{n\A,,A,,...A,,,^^, \n)e-^iE„-,Q„) ^ 0, (5) 

n 

and Furry's theorem will now break down. One may say that the medium, being charged, 
manifestly breaks charge conjugation invariance and these Green's functions are thus 
finite, and will lead to the appearance of new processes in a perturbative expansion. The 
appearance of processes that can be related to symmetry-breaking in a medium has been 
noted elsewhere [Q] . 

Let us, now, focus our attention on the diagrams of Fig. |^ for the case of two gluons 
and a photon attached to a quark loop (the analysis is the same even for QED i.e., for 
three photons connected to an electron loop). Such a process does not exist at zero 
temperature, or even at finite temperature and zero density. At finite density this leads to 
a new source of dilepton or photon production {gg l^l^). In order to obtain the full 



FIGURE 2. The two gluon photon effective vertex as the sum of two diagrams with quark number 
running in opposite directions. 



matrix element of a process containing the above as a sub-diagram one must coherently 
sum contributions from both diagrams which have fermion number running in opposite 
directions. The amplitude for T^pv(= T^P^ + T^P^) is: 



P„„^— (271)-^ {q + k-pYq\q-pY 

Again, the extension of Furry's theorem to finite temperature does not hold at finite 
density: as, if we set n ^ — n — 1, we note that 7^— Qq^ and, as a result, r^P^(/j, T) 7^ 
— r^P'^(/j, r) . Of course. If we now let the chemical potential go to zero ^ 0), we note 
that for the transformation n ^ — n — 1, we obtain qQ —qo, and, thus, r'^P^(0, T) — »• 
— r^P^(0, T) . The analysis for fermion loops with larger number of vertices is essentially 
the same. 



A Realistic Calculation 



To calculate the contribution made by the diagram of Fig. 0, to the dilepton spectrum 
emanating from a quark gluon plasma, we calculate the imaginary part of the photon 
self-energy containing the above diagram as an effective vertex (Fig ^ see reference []T]] 
for details). 



where, the ©'s represent bare gluon propagators, and the T's represent the effective 
vertices. We calculate in the limit of photon three momentum p = 0. The imaginary part 
of the considered self-energy contains various cuts. We concentrate, solely, on the cut 
that represents the process of gluon-gluon to e+e^(see Fig.@). This is a process, which, 
to our knowledge, has not been discussed before. The other possible cuts represent finite 
density contributions to other known processes of dilepton production. 




FIGURE 3. The photon self-energy at three loops and the cut that is evaluated 



The differential production rate for pairs of massless leptons, with total energy E and 
and total momentum p = 0, is given in terms of the discontinuity in the photon self- 
energy M, as 

nw (y I I I 

(8) 



dW ,^ a 1 1 1 „D,^x 



where a is the electromagnetic coupling constant. The rate of production of a hard lepton 
pair with total momentum p = 0, at one-loop order in the photon self-energy (i.e., the 
Born term ), is given for three flavours as 



dEd^p^ ' ^ ' ^ ' ^Tt 

In the above equation, the first term on the r. h. s. is the contribution from the up and 
down quark sector; and the second part is the contribution from the strange sector. 
In a realistic plasma, the net charge and baryon imbalance is caused by the valence 
quarks brought in by the incoming charged baryon-rich nuclei. The baryon and charge 
imbalance is, thus, manifested solely in the up and down quark sector; hence, the 
chemical potential influences only the distribution function of the up and down quarks. 
The strange and anti- strange quarks are produced in equal numbers in the plasma; 
resulting in a vanishing strange quark chemical potential. Thus, dilepton production, 
from the channel indicated by Fig.@, will only receive contributions from the up and 
down quark flavours. 

The initial temperatures of the plasma, formed at RHIC and LHC, have been predicted 
to lie in the range from 300-800 MeV [0, H). For this calculation we use estimates of T = 
400 MeV(Fig. |) and T = 600 MeV(Fig. g). To evaluate the effect of a finite chemical 
potential, we perform the calculation with two extreme values of chemical potential 
/J = O.IT (left plots) and /j = 0.5r(right plots) [Q]. The calculation, is performed for 
three massless flavours of quarks. In this case, the strong coupling constant is (see [|l0|l ) 



The differential rate for the production of dileptons with an invariant mass from 0.5 to 
2.5 GeV is presented. On purpose, we avoid regions where the gluons become very soft. 
In the plots, the dashed line is the rate from tree level qq (Eq. @); the solid line is 
that from the process gg e+e^. We note that in both cases the gluon-gluon process 
dominates at low energy and dies out at higher energy leaving the qq process dominant 
at higher energy. 




FIGURE 4. The differential production rate of back-to-back dileptons from two processes plotted 
against dilepton invariant mass. The dashed line represents the contribution from Born term qq e+e^. 
The solid line corresponds to the process gg e+e^. Temperature is 400 MeV. Quark chemical potential, 
in the first figure, is O.IT The second figure is the same as the first, but, with /j = Q.5T 
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FIGURE 5. Same as above but with T=600MeV. 



High Temperature Limit 



As the reader may have noted, in the above calculation, a ~ 0.3, g ~ 2, thus, we 
are not unequivocally in the perturbative regime. At asymptotically high temperatures 
(T oo), however, ^ ^ 0, in this limit one may make the Hard Thermal Loop (HTL) 
approximation [llTl O] (note: only the main results will be quoted here; the details will 



appear elsewhere D15D). At very high temperature, contributions from loop diagrams are 
dominated by loop momenta (q) of the order of the temperature T. If the momentum 
flowing in the external legs is much smaller (of the order of gT, i.e., q ^ k,p), then loop 



corrections are of the same order in g as the bare diagrams, and, thus, must be resummed 
into the tree amplitudes. 

|-j|2 Ip-op 

In this limit {q + p — k)a — qa, and £'g+p ^Eq + p- q + ^ + . On performing 



the full Matsubara sum and taking the HTL limit for the numerator, we get 

d^q 



-eg^Tr[ffffft] 



qsuaq-s23q-s3,8 



- S3- 



(271)3 2 * II I I I I I 'PO-S2E2+S3E3 
5 1 (n (£3 - 53AJ) - n (£3 + S3IJ) ) + ssjnjEi + sgj) - n{Ei - Sj/j)) 

kO+SiEi+S3E3 

sx{n{E2 - S2H) - n{E2 + S2^i) ) + S2{n{Ei +Sl^l)- n{Ei - si/u) ) 



S2 



k^-pO+SiEi+S2E2 



(11) 



where, = (l,(f 1,^25 <?3): and 4- = i—^Ai^qiAs)- Note that if /j is set to zero this 
amplitude vanishes identically. The conventional HTL term ( i.e., terms proportional to 
(gT)^ ) from triangle graphs such as these would come from the terms with si = +,S2 = 
—,S3 = — and ^'i = —,S2 = +,^3 = +• The term proportional to (gT)^ is 



T, 



gfe^(-fe.4) ^ '-^{{p-k) -q) - '-^{{p-k) -q) 



'-^(-k-q) 



q-,aq-,?,q-,8 

p^+p-q 



k^-k-q 



-k-q) 



k^ + k-q 



+ ■ 



dn{E2+^j) 
dEo 



k^-p^ + ip-k)-q 
{{p-k)-q)-'-^ 



k^-p^-{p-k)-q 



+ 



{{p-k).qy 



Note that if we change q —q in the last two terms then 4+,a ~q-,a and the 
whole term becomes identically zero. Thus there is no term proportional to (gT)^ in this 
diagram. The next term is of order g^T and is nonzero [[T5|]. However, as is well known, 
the HTL approximation of the qqj vertex is of the order of (gT)^ [p^. Thus, in the high 
temperature approximation, the ggy HTL vertex is suppressed as compared to the qqy 
HTL-resummed vertex. As a result the dilepton production rate from resummed two- 
gluon-fusion process will also be suppressed compared to the rate from the resummed 
Bom term, in this limit. 

On summing up all the HTL contributions to the self energy of the gluon we get two 
dispersion relations for the transverse and longitudinal modes of the gluon. A plot of the 
two modes [|T5]] is as shown in Fig. (^. The upper branch is the dispersion relation for 
the transverse quasi-particle excitation, the lower one is the longitudinal excitation. The 
straight line, corresponding to free gluons is shown for reference. Note that there is no 
minima in either branch ( except at = ), unlike in the case of quarks [[l4|]. Hence, we 
will not see any sharp Van Hove peaks [|T^ in the resulting dilepton spectra emanating 
from this process. One may thus expect that the high temperature (and, perhaps, as 
a result, the very low invariant mass) dilepton production spectra emanating from in- 
medium gluon-gluon fusion will be suppressed compared to that from in-medium quark 
anti-quark fusion. However, as we have noted in the previous section the intermediate 




FIGURE 6. HTL resummed dispersion relations of the gluon in a finite temperature, finite density 
medium. 



invariant mass rate from bare gluon-gluon fusion is comparable and may be larger than 
that from bare quark anti-quark fusion. 

The entire calculation above is at full chemical and thermal equilibrium. In the early 
part of an ultrarelativistic heavy ion collision, the gluon number has been predicted to be 
much higher than at full chemical equilibrium. In such a scenario, the contributions to 
dilepton spectrum from processes such as those presented here will probably out-shine 
those from other channels. 
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